Resonance superfluidity in a quantum degenerate Fermi gas 
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We consider the superfluid phase transition that arises when a Feshbach resonance pairing occurs in 
a dilute Fermi gas. We apply our theory to consider a specific resonance in potassium (''"K), and find 
that for achievable experimental conditions, the transition to a superfluid phase is possible at the 
high critical temperature of about 0.5 Tf. Observation of superfluidity in this regime would provide 
the opportunity to experimentally study the crossover from the superfluid phase of weakly-coupled 
fermions to the Bose-Einstein condensation of strongly-bound composite bosons. 



The achievement of Bose-Einstein condensation in 
atomic vapors has given great impetus to efforts to re- 
alize superfluidity in dilute fermionic alkali gases. While 
conditions of quantum degeneracy have been obtained 
in potassium ("^^K) the lowest achievable temper- 
atures to date have been limited to around Q.2Tf [1). 
Although this limit is essentially technical in nature, it 
appears likely that it will be necessary to utilize a strong 
pairing mechanism yielding superfluid transition temper- 
atures close to this value. 

Even in high- Tc superconductors, the typical critical 
temperatures are of the order of fO^^Tp. In the con- 
text of strong-coupling superconductivity there has been 
much work on constructing minimal models to study 
the crossover from the seminal BCS theory Q for con- 
ventional superconductivity to the Bose-Einstein con- 
densation of tightly bound pairs, passing through non- 
perturbative regimes in Tc/Tp [|,|. In this letter, we 
treat explicitly a short range quasibound resonant state 
by extending the theory given in Refs. Q to predict the 
existence of a Feshbach resonance superfluidity in a gas 
of fermionic potassium atoms. This system has an ultra- 
high critical phase transition temperature in close prox- 
imity to the Fermi temperature. This is a novel regime 
for quantum fluids, as illustrated in Fig. |l| where our sys- 
tem and others which exhibit superfluidity or BEC are 
compared. 

The seminal Bardeen-Cooper-Schrieffer (BCS) the- 
ory [Q of superconductivity applied to a dilute gas con- 
siders binary interactions between particles in two distin- 
guishable quantum states, say 1 1) and 1 1). For a uniform 
system, the fermionic field operators may be Fourier- 
expanded in a box with periodic boundary conditions 
giving wavevector-fc dependent creation and annihilation 
operators a^^ and Ofco- for states \a). At low energy, the 
binary scattering processes are assumed to be completely 
characterized by the s-wave scattering length a in terms 
of a contact quasipotential U = Anh'^an/m, where n is 
the number density. The Hamiltonian describing such a 
system is given by 



where = h^k^ /2m is the kinetic energy, m is the mass, 
and the constraint = ki -\- k2 — k^ gives momentum 
conservation. 
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FIG. 1. A log-log plot showing six distinct regimes for 
quantum fluids. The transition temperature Tc is shown as a 
function of the relevant gap energy 2A. Both quantities are 
normalized by an effective Fermi temperature Tp. For the 
BCS systems in region (a), and the systems in the cross-over 
region (b), 2 A is the energy needed to break up a fermion 
pair, and Tp is the Fermi energy. For the systems in region 
(c), which are all strongly bound composite bosons and ex- 
hibit BEC phenomenology, 2A is the smallest energy needed 
to break the composite boson up into two fermions, i.e. ion- 
ization to a charged atomic core and an electron, and Tp is 
the ionic Fermi temperature. 

For a negative scattering length, the thermodynamic 
properties of the gas show a superfluid phase transition at 
a critical temperature Tc which arises due to an instabil- 
ity towards the formation of Cooper-pairs. When the gas 
is dilute, as characterized by the inequality n\a\'^ ^ 1 (or 
equivalently fci?|a| <^ 1 where kp is the Fermi wavenum- 
ber), the application of mean-field theory gives a well- 
known solution for the ratio of Tc to the Fermi tempera- 
ture Tp 0: 
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The exact prefactor to the exponential depends on the 
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precise form of the analytic integral approximations made 
in the derivation. Several papers have pointed out that 
the presence of a scattering resonance in dilute alkali 
gases can be used to obtain a very large negative value 
for the scattering length This promises the oppor- 
tunity for the system to enter the high- Tc superfluidity 
regime as the ratio in Eq. approaches unity. However, 
direct application of the BCS theory close to resonance 
then becomes speculative due to the potential breakdown 
of a number of underlying assumptions: (1) Exactly on 
resonance the theory fails as the scattering length passes 
through ±00 and the Hamiltonian in Eq. ([l]) can not be 
defined, (2) For the mean field approach to be accurate it 
is required that there be many particles inside a volume 
associated with the spatial scale of a Cooper-pair. This 
condition begins to break down as Tc approaches Tp, 
(3) The theory of the dilute gas is formulated on a per- 
turbation approach based on an expansion in the small 
parameter njap. When this parameter approaches unity 
the perturbation theory fails to converge. These points 
show that care should be taken in applying Eq. (g) near 
the point of resonance where the basis for the conven- 
tional mean-field theory is not well-founded. 

Despite these limitations, on general grounds, one 
would expect to be able to derive a renormalizable low- 
energy effective field theory even in close proximity to a 
resonance. This statement is based on the identification 
that at relevant densities the range of the interparticle 
potential is always orders of magnitude smaller than the 
interparticle spacing. Here we present a theory of super- 
fluidity in a gas of dilute fermionic atoms which handles 
correctly the scattering resonance and places the tran- 
sition temperature to the superfluid state in the experi- 
mentally accessible range. 

While the scattering length a usually characterizes the 
range of the interatomic potential for a collision, this is 
a poor approximation in the vicinity of a scattering res- 
onance. The scattering properties are completely deter- 
mined by the positions of the bound states in the interac- 
tion potentials. In a multichannel system, a bound state 
may cross the threshold as a function of magnetic field 
and enter the continuum, resulting in a field-dependent 
Feshbach scattering resonance As this occurs, the 
scattering length becomes strongly dependent on the 
field, and exactly at threshold it changes sign by pass- 
ing through ±(X). 

When such resonance processes occur, it is necessary 
to formulate the Hamiltonian by separating out the res- 
onance state and treating it explicitly. This is motivated 
by the microscopic identification of two types of scatter- 
ing contributions: one from the scattering resonance, and 
one from the background non-resonant processes that in- 
cludes the contributions from all the other bound states. 
The non-resonant contributions give rise to a background 
scattering length Obg which is a good characterization of 
the potential range. The corresponding quasipotential 
in that case is given by C/bg — ■iTrh^a^jgn/m. The Fesh- 
bach resonance occurs due to a coupling with a molecular 



state, that is long-lived in comparison with characteristic 
non-resonant collision timescales. This state is a com- 
posite boson which is described by bosonic annihilation 
operators 6^,. It is parameterized by a detuning energy 
from threshold, denoted by 2iy, that is dependent on the 
value of the magnetic field. The coupling strength of b,^ 
to the two-particle continuum is well characterized by a 
single coupling constant g, independent of k. These con- 
siderations imply that the Hamiltonian given in Eq. (Q) 
is not sufficient to account for the important resonance 
processes and must be extended to incorporate explicitly 
the coupling between the atomic and molecular gases: 

fe k 

-I- C/bg "LT"Li«fe3i"fc4T 

ki ...k^ 
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Evolution generated by this Hamiltonian conserves the 

particle number N = Efc(alT"fcT +«ili«fei) +2Efc/'fc^fc- 
Note that the Hamiltonian does not contain a explicitly, 
and that the field dependence of the scattering is com- 
pletely characterized by the parameters: g, u, and t/bg- 
The magnitude of g is derived in the following way. We 
define k as the product of the magnetic field width of the 
resonance and the magnetic moment difference of the Fes- 
hbach state and the continuum state. For large values of 
h', the boson field 6^ can be adiabatically eliminated from 
the theory, and then g — nUhg is required in order for 
the scattering properties to have the correct dependence 
on magnetic field JlOt . 

The essential point is that this Hamiltonian, founded 
on the microscopic basis of resonance scattering, is well- 
behaved at all detunings v; even for the pathological case 
of exact resonance. The diluteness criterion is now given 
by constraints which require both the potential range 
and the spatial extent of the Feshbach resonance state, 
to be much smaller than the interparticle spacing (e.g. 
n|abgP < 1). 

We apply this Hamiltonian to derive the self-consistent 
mean-fields for given thermodynamic constraints by for- 
mulating a Hartree-Fock-Bogoliubov theory The 
mean-fields present include the fermion number / = 

X]fc(«fct'^fc|)j tl^e molecule field </),„ = {b^^^) taken 
to be a classical field, and the pairing-field p — 



) ||13]. It is well-known that such a theory 



^ki^k^'^-kl ^ 

must be renormalized in order to remove the ultraviolet 
divergence which arises from the incorporation of second- 
order vacuum contributions. This implies replacing the 
physical parameters in the Hamiltonian, U, g, and by 
renormalized values so that observables are independent 
of a high momentum cut-off used in the formulation of 
the effective field-theory |13 . In order to diagonalize the 



Hamiltonian, we construct Bogoliubov quasiparticles ac- 
cording to the general canonical transformation [n4| 
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Given single particle energies, Uk — Ck — fJ^ + Uf, where 
fj, is the chemical potential, and the gap parameter in 
the quasiparticle spectrum A — Up — g4>m, the two 
transformation angles are specified as tan(26') = |A|/C//j 
and = \4>m \ exp(?7). The corresponding quasiparticle 
spectrum is Ek — \/U'^ + A'^. Dropping terms of higher 
order than quadratic in the fermion operators, gives the 
resulting many-body Hamiltonian 

H - fiN = 2(1^ - ^i)\cj)„A^ 

+ E(^fc+^^-KT"feT+4i«fei-l))- (5) 
fe 

which is now in diagonal form. 

The next task is to calculate the thermodynamic so- 
lutions. Equilibrium populations for the quasiparticles 
are given by the Fermi-Dirac distribution. The fermion 
number and pairing field are not only inputs to the 
Hamiltonian, but also determine the quasiparticle spec- 
trum. Therefore, they must be self-consistent with the 
values derived by summing the relevant equilibrium den- 
sity matrix elements over all wave numbers. In prac- 
tice, at a given temperature, chemical potential, and 
molecule number 0™, this requires an iterative method 
to locate self-consistent values for / and p. The value 
of 0m is calculated by minimizing the grand potential 
$(3 = — fc^TlnS at fixed temperature and chemical po- 
tential, with kb denoting Boltzmann's constant. The par- 
tition function S = Tr[exp(— (iJ — fj,N)/ki,T)] is found 
from Eq. (|^). This procedure is mathematically equiv- 
alent to minimizing the Helmholtz free-energy at fixed 
temperature and density and corresponds uniquely to the 
maximum entropy solution. This solution has an associ- 
ated particle number, (iV) = — 9$G/<9/i taken at con- 
stant temperature and volume, which must match the 
actual particle density of the gas, so that the final step 
is to adjust the chemical potential until this condition is 
satisfied. The whole procedure is repeated over a range of 
temperatures to determine the locus of thermodynamic 
equilibrium points. For large positive detunings, where 
the molecule field could be eliminated from the theory en- 
tirely, regular BCS theory emerges. For this case, when 
the scattering length a is negative the behaviour of the 
critical temperature on 1/a is given by the usual expo- 
nential law . 

As an example of the application of this theory, we 
study the experimentally relevant system of fcrmionic 
*'^K atoms equally distributed between the two hyper- 
fine states which have the lowest internal energy in the 
presence of a magnetic field. The values of our interaction 
parameters abg = 176 ao and n/ki, = 657 /xK are obtained 
from jl^. We fix the total density to be n = 10^^ cm""^, 
a typical experimental value expected for this quantum 
degenerate gas in an optical trap. We set the detuning 
to be = +Ef so that the quasi-bound state is detuned 



slightly above the atomic resonance. For a temperature 
above T^, the grand potential surface is shaped like a 
bowl, and the value of (pm which minimizes the grand 
potential is (pm = 0, associated with the self-consistent 
solution p = 0. For T < Tc, the grand potential surface 
is shaped like a Mexican hat, and its minimum is given 
by a 0™ with non-zero amplitude and an undetermined 
phase. The superfiuid phase transition therefore leads to 
a spontaneously broken symmetry. The value of can 
be clearly found from Figs. || and ||, where we show the 
chemical potential, the molecular density, and the gap 
as a function of temperature. We find for our parame- 
ter set for "'"K and almost zero detuning a remarkably 
high value for the critical temperature Tc ~ 0.5 Tp, i.e. 
Tc ~ 0.6 fiK. Furthermore we find a weak dependence of 
Tc « 0.5 Tp on the density, so that the value of Tc has 
more-or-less the same density behavior as Tp. When we 
increase the detuning to v = +17.6 Ep (this corresponds 
to a magnetic field detuning of 0.5 G away from the Fcs- 
hbach resonance), the value of Tc drops to approximately 
0.25 Tp. 
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FIG. 2. Chemical potential as a function of temperature 
for the system of resonance pairing (solid line). The second 
order phase transition occurs at Tc ~ 0.5 Tf where a clear 
cusp is visible. The dashed line shows the chemical potential 
of a non-interacting Fermi gas. 

The system of '*°K atoms, equally distributed among 
the two lowest hyperfine states, is a good candidate for 
demonstrating the superfiuid phase transition. It not 
only exhibits a Feshbach resonance, but also, the in- 
elastic binary collision events are energetically forbidden. 
Three-body interactions are highly suppressed, since the 
asymptotic three-body wave function should consist of a 
product of three s-wave two-body scattering wave func- 
tions. In a three-body interaction, two-particles are al- 
ways in the same initial hyperfine state, and therefore the 
corresponding s-wave state is forbidden. The only three 
body relaxation could come from asymptotic p-waves, 
but these have very little contribution at the low tem- 
peratures considered. Although the detailed three-body 
collision problem is an intricate one, this asymptotic sta- 
tistical effect should lead to a large suppression of the 
vibrational relaxation of quasi-bound molecules. 
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FIG. 3. The temperature at the phase transition is also 
visible from the amplitude of the molecular field. This am- 
plitude is non-zero only when the broken symmetry exists in 
the region T < Tc. For T = 0, the molecules form a Bose 
condensed fraction of 1.5% of the total gas sample. The inset 
shows the behaviour of the gap A = Up — g<j>m- The criti- 
cal temperature Tc can be related to the value of the gap at 
T = 0. For comparison, in superconductors the analagous 
gap is simply the binding energy of a fermion pair. 

Current experimental techniques for ultracold gases do 
not produce samples which are spatially uniform. An 
optical dipole trap may be needed to confine the high 
field seeking atoms, and the conditions for the superfluid 
phase transition would be satisfied first in the trap cen- 
ter where the density is highest. The presence of the 
quasi-bound molecules may be a very useful aspect al- 
lowing direct observation of the phase transition through 
imaging the molecular field. 

In conclusion, we have shown that resonance pairing 
in an alkali gas yields a quantum fluid that can undergo 
a superfluid phase transition at a temperature compara- 
ble to the Fermi temperature. This extraordinary prop- 
erty places this system in a regime which lies in between 
BCS-like superconductors, and bosonic systems which 
may undergo BEC. Since the transition temperature is 
larger than the lowest temperatures already achieved in 
a degenerate Fermi gas, it should be possible to study 
this new type of quantum matter and to quantitatively 
compare with our predictions. 
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Note Added: A similar treatment has recently been 
proposed by Timmermans et al. |l6|. 
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